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Abstract 



In this article, we propose a resolution to the paradox of apparent superluminal velocities 
for tunneling particles, by a careful treatment of temporal observables in quantum theory 
and through a precise application of the duality between particles and waves. To this end, 
Q_|. we employ a new method for constructing probabilities associated to quantum time measure- 

ments that provides an explicit link between the tunneling time of particles and the associated 
quantum fields. We demonstrate that the idea of faster-than-light speeds in tunneling follows 
from an inadmissible use of classical reasoning in the description of quantum systems. Our 
results suggest that direct measurements of the transit time in tunneling could provide a new 
testing ground for the predictions of quantum theory versus local hidden-variables theories. 

> 

1 Introduction 

in 

One of the most fundamental features of quantum theory is the duality between particles and waves, 
namely, the fact that the same quantum system can exhibit particle-like or wave-like properties 
in different experimental configurations. The particle and the wave aspects of any system are 
complementary, they cannot be combined into a single unified description of microscopic properties 
that is independent of the experimental set-up. In this article, we show that the long-standing 
^ | paradox of apparent super-luminal velocities in quantum tunneling can be resolved by a careful 
implementation of particle- wave duality in the context of Quantum Field Theory (QFT). 

The superluminality paradox in tunneling originates from the early days of quantum theory [1] , 
in the effort to identify how long it takes a quantum particle to tunnel through a potential barrier. 
The search for an answer to this question has led to a large number of candidates for the tunneling 
time [2j 13] rather than to a single expression derived unambiguously from first principles. Many 
existing definitions imply that tunneling times saturate in the opaque-barrier limit, thus suggesting 
superluminal speeds for particles traversing sufficiently long barriers (the Hartmann effect [1]). 

Several experimenters have reported superluminal tunneling velocities in electromagnetic ana- 
logues of quantum tunneling [5J El El El El EE]- The analogy is based on the mathematical corre- 
spondence between the classical Helmholtz equation for the electromagnetic field in inhomogeneous 
dielectric media and the time-independent Schrodinger equation in presence of a potential. The 
experiments above measure the group delay t d associated to electromagnetic pulses crossing a bar- 
rier of length d. In classical field theory, the group delay is standardly interpreted as the transit 
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time of a signal through a medium or a device. Hence, the group velocity d/t^ (which is found 
greater than the speed of light) is interpreted as the transit velocity for waves crossing the barrier. 
However, no direct measurement of the transit time is involved in these experiments; the transit 
time is inferred from the group delay. 

Our approach to the superluminality paradox is based on the observation that the relation 
between group delay and transit time is a feature of classical physics and not of quantum physics. 
In quantum theory, there is no such relation, because group delay and transit time are incompatible 
observables that are determined in distinct types of experiment. The former is an observable in 
experiments that measure the wave aspects of a quantum system, while the latter refers to the 
propagation of localized excitations of quantum fields, i.e., particles. The interpretation of group 
delay as a transit time violates a fundamental rule of quantum mechanics, Bohr's complementarity 
principle, according to which "evidence obtained under different experimental conditions cannot 
be comprehended within a single picture" [TT] . 

Group delay in electromagnetic analogues of tunneling has a natural interpretation as the 
lifetime of energy stored in the barrier [3j [12]. To address the issue of superluminality, we must 
consider experiments in which the particles' transit times (or times of arrival) are direct physical 
observables. Such experiments are well understood in classical theory. They typically involve 
a particle source and a particle detector separated by distance L at rest with respect to each 
other. Source and detector are equipped with a pair of synchronized clocks. The transit time is 
defined as the difference t between the clock readings of detection and emission respectively, and 
the traversal velocity is defined as L/t. When considering quantum particles, the time difference 
t becomes a random variable that may take different values in different runs of the experiment. 
Thus, a quantum description requires the construction of a probability density P(L, t) for the 
transit time. When a potential barrier is placed along the line connecting emitter and detector, all 
information about the temporal aspects of tunneling is contained in the probability density P(L, t) 

m. 

However, the construction of a probability density P(L,t) for transit times associated to tun- 
neling experiments is a highly non-trivial task. The derivation of probabilities with respect to time 
in quantum theory has run into the same type of problems as the identification of tunneling time 
|14j . These problems originate from the subtle role of time in quantum theory, where the time t 
appearing in Schrodinger's equation is an external parameter and not an ordinary observable, like 
position or momentum. This implies that the squared modulus of the time-evolved wave-function 
\ip(x,t)\ 2 is not a density with respect to t, and, hence, it cannot serve as a definition for the 
required probabilities. 

A significant advance towards the resolution of the issues above is the development of an 
algorithmic method for constructing Quantum Temporal Probabilities (QTP) associated to any 
experimental configuration [TJ] — see also Ref. [18]. The QTP method incorporates the detector 
degrees of freedom into the quantum description, so that the temporal probabilities are always 
defined with respect to specific experimental set-ups. The key idea is to distinguish between the 
roles of time as a parameter to Schrodinger's equation and as a label of the causal ordering of 
events [16] . This important distinction allows for the definition of quantum temporal observables. 
In particular, we identify the time of a detection event as a coarse-grained quasi-classical variable 
[TT] associated to macroscopic records of observation. Thus, the time variables in QTP correspond 
to macroscopic observable magnitudes, such as the coincidence of a detector 'click' with the reading 
of a clock external to the system. 

In Ref. [19], we used the QTP method in order to derive of a probability density P(L,t) for 
the time t of particle detection in tunneling set-ups. The derivation involves a QFT description of 
the tunneling process that is necessary for a definitive resolution of the superluminality paradox. 
The reason is that superluminal speeds imply a violation of the principle of local causality, and 
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this principle is implemented in quantum theory only if the interactions are expressed in terms of 
local quantum fields [20] . 

In what follows, we briefly present the relevant results of Ref. [19]) (Sec. 2), and then we use 
these results in order to propose a resolution of the superluminality paradox (Sec. 3). 



2 The QTP method 

In Ref. [19], we constructed the detection probability for quantum particles of mass m and charge 
e , described in terms of a quantum field <fi(x) tunneling through a potential barrier. Considering 
for simplicity a single spatial dimension, we assume that the barrier is localized in a region D = 
[—d/2, d/2], where d is the length of the barrier. A particle detector is localized around x = L » d. 
In the QTP method, the explicit quantum modeling of the particle detector leads to a unique 
probability density P(t) associated to an experiment. The most important element in the modeling 
of the detector is the specification of the interaction Hamiltonian Hj between the detector and the 
quantum field 4>(x). The requirements of Lorentz covariance and unitarity imply that Hj must be 
a local functional of the fields <f>(x) and <jy{x) [2D]- We select 



H I= dx 



(j)(x)®J j (x) + ft(x)®J(x) , (1) 

where J and are current operators on the Hilbert %det of the detector degrees of freedom. The 
initial state of the detector |^ ) is assumed to satisfy the condition J(x)\^q) = 0, which guarantees 
that the detector is sensitive to particles rather than anti-particles. 

For the set-up described above, we derived a general expression for the probability density 
P{L,t) 

P(L,t) = J dr J dzR(z,T){i> \ft(L-^t-^(L+^t+^)\i; ). (2) 

where l^o) is the initial state of the field and H the corresponding Hamiltonian operator. All 
information about the detector is contained in a kernel R(z, t) that smears the field correlation 
function. It is important to emphasize that Eq. (j2J) requires no specification of the Hamiltonian H 
or the quantum field 0. In particular, it does not depend on how we choose to model the potential 
barrier. It is valid for a barrier modeled by a background classical field, but also for barriers 
defined in terms of the full quantum interaction of the field </> with other fields (i.e., the quantum 
electromagnetic field). The only assumption in the derivation of Eq. (|2j) is that the detector at 
x = L detects particles corresponding to the field 0. Thus, Eq. (j2j) reveals a relation between 
probabilities for detection time and field correlation functions that persists in any quantum field 
theory. 

For concreteness, we specialized to the case that the potential barrier corresponds to a back- 
ground static electromagnetic potential A^x) = (A (x),0), where A (x) differs from zero only in 
the spatial region D. Then, the classical Hamiltonian for the field <fi{x) is 



H 



J dx [\tt\ 2 + \d x <f)\ 2 + m 2 \<P\ 2 + tV{rr<j) - tt*0*)] , (3) 



where ir(x) is the conjugate momentum to 4>(x) and V(x) = eA (x). The quantized Hamiltonian H 
is expressed in terms of the standard creation and annihilation operators: a(x), at(x) for particles 
and b(x),b^(x) for anti-particles 
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H = j dxdx' ya jf (x)hi(x, x')a(x) + P(x)h,2(x, x')b(x)J , (4) 

where h\{x,x') and h,2(x,x') denote matrix elements of Hamiltonian operators hi and h 2 on the 
Hilbert spaces Hi and H* associated to single particle and a single anti-particle, respectively. The 
operators hi and hi are defined as hip — ho ± V( x ), where ho = a/— d 2 + w^ 2 is the Hamiltonian 

for a sine le free particle, and V := ± (^V/iq 1/2 + ^V/i" 1/2 ) = ^ + [[V, /iq /2 ], /iq 1/2 ] is the 
QFT-corrected potential term. 

For a positive- valued potential V(x), the single-particle Hamiltonian /ii has continuous spec- 
trum. There exists a pair of orthogonal energy eigenstates to each value of energy E > m. We 
denote such eigenstates as f k+ and f k -, where k = y/E 2 — m 2 > 0. The eigenstates f k+ correspond 
to only positive momentum flux at x = oo. They are of the form 



1 < p ikx I U, p -ikx ^ _d 

M } ~ V2^ I T k e ikx x>l [b) 

where T k and R k are the usual transmission and reflection coefficients for a right-moving plane 
wave of momentum k. We parameterize the eigenstates f k -( 

fk-{x) = 1 \-w k f k+ (x) + f k +{-x)] , (6) 



w l 



where w k = \{T k R* k + T k Rk) is essentially the overlap between f k+ and its parity transform. For 
parity-symmetric potentials, V(x) = V(—x), the coefficient w k vanishes. The Hamiltonian h% for 
a single antiparticle, has discrete-spectrum eigenstates g n {x) for energies E n < m, in addition to 
the continuous spectrum eigenstates g k (analogous to the f k above) for E > m. 

Defining a k ± = J dxa(x)f k± (x), b k ± = J dxb(x)gl ± (x), and b n = J dxb(x)g^(x), we express the 
Heisenberg-picture field operators as (j>(x,t) = <fii(x,t) + (f>2(x,t), where 



*) = E / -4=^hA*)e- iEkt (7) 

Jo V ^E k 

+ E^^y^ (8) 



Jilar. 



We assume that the particles are emitted from a source, localized around x = —x < —d/2, 
in an initial state l^o) — / dxa) (x)ipo{x)\0) , where if>o(x) is a single-particle wave-function with 
positive momentum and |0) is the field vacuum. Substituting Eqs. ([TUB]) into Eq. ([2D, we find that 
the contribution from the anti-particle component 02 is negligible for L » d. For the detector 
far from the barrier, the probability density P(L < t) becomes (15] 



2 

(9) 



where v k = kj\/k 2 + m 2 is the relativistic velocity and tpo(k) is the initial wave-function in the 
momentum representation. In Eq. fl9]), all information about the detector is encoded in the 



P(L,t) 



— a(k)\v k \A k tpo{k)e 

Z7T 



ikL—iEkt 
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absorption coefficient a(k) for particles of momentum k, and all information about the barrier is 
contained in the complex amplitude 



A k 



T k — w k R k 



(10) 




In absence of the potential barrier, A k = 1. Then, Eq. reduces to the time-of-transit 
distribution for free relativistic particles derived in Ref. [T5] . 

3 Eliminating the superluminality paradox 

We consider the probability distribution P(L, t), Eq. OH]) for an initial state tp localized at X — — Xq 
and with a narrow momentum distribution centered around k. In absence of the barrier, P(L,t) 
is sharply localized around the mean value (L + xq) / v k of the transit time. The insertion of the 
barrier induces a delay td{k) to the mean transit time, where 



For parity-symmetric potentials, A k = T k . The time delay of Eq. (ITT]) then reduces to the 
phase delay time of Bohm and Wigner [211 122] . Note that Eq. (fTTj) follows from a stationary phase 
approximation to the probability density P(L, t) of Eq. (J9J) that is valid only if P(t) has a single 
maximum. 

In general, the temporal properties of the tunneling process are encoded in the full probability 
distribution P(L,t) associated to a particular experiment, and not in a single parameter like the 
delay time Eq. (ITTj) . The structure of P(L,t) is the only criterion for whether the delay time is 
meaningfully defined, or for the existence of other, physically significant, timescales. Indeed, if the 
probability distribution P(L,t) is not characterized by a single peak, the identification of td{k) as 
a delay time is misleading, if not downright wrong [T9] . 

Nonetheless, Eq. (iTTj) applies to several important cases, such as tunneling through a square 
barrier. In this case, the potential is parity symmetric. This implies that A k = T k , and the time 
delay Eq. ( ITTj) reduces to the standard phase delay time [2J . Our results then still face the challenge 
of potential superluminal signals. For an opaque square barrier, we recover the standard form for 
the delay time tdik) = —d/v k + F{k) [2], where F(k) is a function only of the incoming particle's 
momentum k. If one defines the time r that the particle spent inside the barrier as r = + d/v k , 
then t is independent of the barrier length d. This implies superluminal traversal velocity for long 
barriers. 

The problem with this line of reasoning is that the time delay td{k) is not a genuine quan- 
tum observable. It does not correspond to a self-adjoint operator (or a positive-operator-valued 
measure), and it does not correspond a measurement record in the experiment. The measurement 
records correspond to the time of transit t, not to td- The time delay is a parameter of the proba- 
bility distribution for the detection times and it can only be identified after the probability P(L, t) 
has been determined. In fact, the identification of the the time delay td(k) requires the combina- 
tion of data from two different time-of-transit experiments, one with and one without the barrier. 
Each experiment records a probability distribution for the transit times and td(k) is identified as 
the difference of the corresponding mean values. For this reason, an inference of superluminality 
from the values of the delay time td(k) involves an inadmissible treatment of measurement records 
that violates the complementarity principle. In no experiment is an actual superluminal signal 



t d (k) 



— Im 

v k 



dlog A k 
dk 



(11) 
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recorded. The events of particle emission and particle detection are time-like separated. Their 
proper distance As 2 = [(L + x )/vk + t d ] 2 — (L + x ) 2 is positive in the regime L » d where Eq. 
(TlTI) applies. This is because for L » d the particle spends considerably more time outside the 
barrier region. 

The positivity of As 2 for L » d relies on the fact that for the square barrier, 

t d (k) > -d/v k . (12) 

However, the conclusion is not restricted to this particular case. For L » d, tunneling can be 
described as a scattering process, and the derivatives of the scattering phase shifts must have a 
lower bound in order for the S'-matrix to be compatible with causality [22j [23] . The lower limit in 
the phase shifts implies a lower bound to the delay time t d (k). For non-relativistic particles and in 
absence of bound states, the lower bound Eq. (fT21 applies [21]. In general, the corrections to this 
lower bound are of the order of the particle's de-Broglie wave-length divided by v p , i.e., they are 
microscopic. They do not affect the argument about the time-like separation between the events 
of particle emmission and particle detection. 

Thus, the only conceivable way of recording a superluminal signal would be in an experimental 
set-up where both emitter and detector are placed very close to the barrier so that As 2 < 0. 
This implies that r/v p + td{p) < r, where we write r = L + x f° r the source-detector distance. 
Since ta(p) > —d/v p , the necessary condition for As 2 < is that d < r < r/(l — v p ). The 

velocity v p is bounded above by y^l — (1 + j%)~ 2 '- for larger velocities there is no tunneling. The 
approximation of a background electric field fails for Vo > m, because in this case the Hamiltonian 
Eq. (@J would possess negative-energy eigenstates. This implies an upper bound to the velocity 
v p < y/3/2 ~ 0.87. A less stringent upper bound to the velocity (v p < 2\/2/3 ~ 0.94) follows from 
the requirement that the particle's kinetic energy is less than 2m, so that particle-antiparticle pairs 
are not created spontaneously. Thus, we conclude that the only conceivable range of values for r 
where a superluminal might be possible on the basis of Eq. f fTTj) is 

d < r < 7.5d, (13) 

i.e., the emitter-detector distance r is of the order of the barrier length d. However, in this case 
the stationary phase approximation that leads to the delay time Eq. ( II ip is invalid. The position 
spread a x of the initial wave-packet must be much smaller than r. To see this, note that the spread 
<j x in position of the initial wave-packet must be much smaller than r. For r ~ d, this implies 
that <j x » d. Hence, the momentum spread a p must be very large, a p » 1/d, a condition that 
is incompatible with the stationary phase approximation employed in the derivation of Eq. (jTTj) . 
The probability distribution P(L, t) is strongly deformed and there is no meaningful definition of 
a delay time tdik). 

In fact, the absence of superluminal signals is guaranteed by Eq. ([2]) that relates the probability 
distribution P(L,t), to the correlation functions of a local quantum field. A single-particle state 
\ipo) can be expressed as \ipo) = J dx(f>i(x) f* (x)\0) , for some localized function f(x). For a single- 
particle initial state localized around x = —xq < —d/2, Eq. (T2]) implies that the probability P(L, t) 
is proportional to the four-point function 

(0|<M-x , 0)#(L, t)ML, t)0\(-x o , 0)|0), (14) 

modulo spatial smearing at the points of emission and detection. Since the points x = —Xq and 
x = L are outside the barrier region, the field 0i is effectively free at these points, and they relate 
to the asymptotic fields defined at t — > ±oo via evolution through the free field Hamiltonian. 
In particular, 4>\(— x ,0) corresponds to the in-field and (p 2 (L,t) to the out field of the S'-matrix 
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formalism. For free fields, the four point function Eq. ( 1T%|) factorizes into a product of two point 
functions. Hence, 

P(L,t)~|A(-a;o,L;t)| 2 (15) 

where A(x,x r ;t) = [4>i(L, t), 4>\(— xo, 0)] = (0|</>i(L, t)<j)\( — x , 0)|0). This means that the probabil- 
ity P(L,t) involves propagation of the initial state through the two-point function A(x,x';t) of 
the quantum field. In any causal QFT, this Green's function, constructed from the in and out 
fields must vanish outside the light- cone, for otherwise it would lead to a non-causal S-matrix. It 
follows that the detection probabilities cannot involve superluminal signals. Thus, we claim that 
the QFT description of tunneling time, enabled through the use of the QTP method, eliminates 
the superluminality paradox. 

The same point of principle holds for the electromagnetic analogues of tunneling, mentioned 
earlier, where superluminal group velocities have been recorded. The analogy of those experiments 
to quantum tunneling is based on the following correspondence. The classical Helmholtz equation 
for an electromagnetic field mode at frequency u in an inhomogeneous dielectric medium of 
refraction index n(x) is 

d%E u + {n{x)ujf E^ = 0. (16) 
This is formally analogous to the time-independent Schrodinger equation in presence of a potential 

d 2 x i) + [2m(E - V(x))} = 0, (17) 

1 /2 

if we set n(x)cj = [2m[E — V^(x))] 1 . Depending on the dielectric, it is possible to have evanescent 
waves, which decay very much like the quantum mechanical wave functions in tunneling. 

However, the analogy between the Helmholz and the Schrodinger equations holds only if they 
are viewed as describing classical waves equations. At the quantum level there is a significant 
difference. Time evolution according to Schrodinger's equation is always unitary. However, the 
effective Hamiltonian operator for the quantum electromagnetic field in presence of a dielectric is 

H e ff = ^2n(uj a )uj a a ] a a a) (18) 

a 

where a, a) are creation and annihilation operators of the electromagnetic field, and a labels the field 
modes. Evanescent waves appear for imaginary values of the refraction index n(u a ), in which case 
the operator (1181) becomes non-hermitian. The effective quantum evolution according to Eq. (JTHJ) 
would then be non-unitary. This is to be expected, since the microscopic mechanism generating an 
imaginary refraction index is the absorption of photons by the atoms of the medium. This implies 
that the electromagnetic analogues of tunneling are in fact analogues of time-of-arrival experiments 
in open quantum systems, where the effects of dissipation and noise have to be incorporated in 
the quantum description. 

A full treatment of this effect requires the application of the QTP method for the treatment of 
temporal measurements in open quantum systems, which will be the topic of a different publication. 
Here, we point out that a fully quantum treatment of the electromagnetic field in inhomogeneous 
media requires the complete Quantum Electrodynamics (QED) Hamiltonian for the interaction 
between the electromagnetic field and the medium. That is, we consider a Hilbert space H to t = 
Hem ® %med ® Happ, where Hem is associated to the electromagnetic field, H m ed to the particles 
forming the dielectric medium and T-L ap p associated to the measurement apparatus. The dynamics 
on Hem <8> Hmed is governed by the QED Hamiltonian 

Hqed = H EM ® 1 + 1 ® H med + / dxA\x) ® Ux), (19) 
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where Hem is the Hamiltonian for the free electromagnetic field, H me( i is the self-Hamiltonian for 
the particles in the medium, A l (x) the electromagnetic potential and ji the current associated to 
the medium. 

The coupling of the electromagnetic field with the detector is governed by an interaction Hamil- 
tonian Hj analogous to Eq. (jSJ) 

Ei = J dxA l (x)Ji{x), (20) 

where Ji(x) is the electric current associated to the detector degrees of freedom. We can then 
derive an equation for the probability density of detection P(L, t) analogous to Eq. (j2J) 

P(L,t) = JdrJ dzR t \z,r)^ ,n\Al(L-^t-^)A 1 (L + ^t+^ ,n), (21) 

where R^ij, z) is a kernel incorporating the detector degrees of freedom, and j^o, ^) = IV'q) ® 
where \ipo) is the initially prepared electromagnetic field state, and \Q the initial state of the 
particles in the medium (for example, an energy eigenstate). Using the same arguments leading to 
the estimation Eq. (1151) for P(L, t), we can show that the propagation of the electromagnetic field 
signal in Eq. (|2"T|) is guided by the two-point function of the field A, that corresponds to photon 
emission and detection. If QED is a consistent quantum field theory these functions must be causal, 
irrespective of the initial state |0) of the particles in the dielectric medium. Thus, no superluminal 
signal is to be expected in transit-time measurements of photons through an absorbing medium. 



4 Conclusions 

To summarize, our proposed resolution to the superluminality paradox is the following. The time 
delay td{k) can indeed be determined in some time-of-transit experiments on tunneling set-ups. 
However, it is not a genuine quantum observable, but a parameter of the probability distribution 
P(L,t). An inference of superluminal velocities in the tunneling region (or indeed any velocity) 
from the value of td{k) involves classical reasoning that is incompatible with the rules of quantum 
theory. All physical information about signal propagation is contained in the probability distri- 
bution P(L,t), Eq. (j2J), which is guaranteed to be causal, because, the QTP method enables its 
definition in terms of the correlation functions of a local quantum field theory. 

In this sense, the apparent superluminality paradox in tunneling is a non-trivial manifestation 
of the particle-wave duality in quantum theory. Group velocities may be greater than the speed of 
light, but they are defined in terms of set-ups that measure wave properties of the quantum field, 
and not transit times. The latter correspond to the particle aspects of the quantum field, and they 
are measured in different experiments. Inferences of superluminality follow from attempts to relate 
observables defined in different experiments. Indeed, the only way to relate physical magnitudes 
associated to different experiments is through the introduction of hidden variables, i.e., variables 
describing microscopic properties of the system that are not contained in the formalism of standard 
quantum theory. But superluminality would hardly be surprising in this case. Superluminal 
velocities are to be expected in any hidden-variables theory that reproduces the predictions of 
quantum mechanics, by virtue of Bell's theorem [25J. 

For this reason, we believe that time-of-transit experiments in tunneling systems could be 
of great significance for the foundations of quantum mechanics. We expect that the probability 
distribution Eq. (Q for the time of transit is incompatible with any local hidden-variable theory, 
i.e., a hidden- variable theory that admits no superluminal propagation. Thus, experiments aiming 
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to the confirmation of Eq. could provide a new ground for testing the predictions of quantum 
theory Vs. local realism. They could demonstrate that quantum 'non-locality' refers not only to 
the correlations between quantum subsystems, but also to the values of temporal observables in 
non-composite quantum systems. 
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